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Two of the most commonly used tools to constrain the primordial non-Gaussianity are the bis- 
pectrum and the Minkowski functionals of CMB temperature anisotropics. These two measures of 
non-Gaussianity in principle provide distinct (though correlated) information, but in the past con- 
straints from them have only been loosely compared and not statistically combined. In this work 
we evaluate, for the first time, the covariance matrix between the local non-Gaussianity coefficient 
/nl estimated through the bispectrum and Minkowski functionals. We find that the estimators 
are positively correlated, with correlation coefficient r ~ 0.3. Using the WMAP7 data to combine 
the two measures and accounting for the point-source systematics, we find the combined constraint 
/nl = 12 ± 29. 



Introduction. Detection of any departures from Gaus- 
sianity in the distribution of primordial fluctuations 
would give important information about inflation. Pri- 
mordial non-Gaussianity (henceforth NG) imprints sig- 
natures on the cosmic microwave background (CMB) and 
large-scale structure, and these cosmological probes can 
in turn provide excellent constraints on primordial NG 
and thus inflationary models; for reviews, see [IH3]. 

Two of the principal statistics on the CMB used to 
constrain NG are the bispectrum (harmonic transform 
of the three-point correlation function) of the CMB tem- 
perature fluctuations, and Minkowski functionals (hence- 
forth MF) which roughly measure the connectedness or 
morphology of the CMB field. In the "local" model of 
NG, the primordial curvature perturbation $(x) has a 
quadratic term correction: $ = $ G + /nl($g — (^g))> 
where <&g is an auxiliary Gaussian field [5]. Recent con- 
straints obtained on the non-linear coupling constant /nl 
using the WMAP data are /nl = 37 ± 20 from the bis- 
pectrum analysis [6] (see also PJE]), an d /nl = 20 ± 42 
from the MF analysis [S]- 

Since MF are morphological statistics, they probe NG 
both in configuration space and to all orders of the statis- 
tics of the temperature anisotropy field. This means 
that they sample the anisotropy map differently from the 
usual bispectrum (and higher order polyspectra) mea- 
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surements, albeit in a suboptimal way - this fact is cru- 
cial as joint constraints will in principle yield different 
constraints. Furthermore, unlike the bispectrum estima- 
tors which require a template (i.e. fc-space configuration 
with a free amplitude) such as the local, equilateral or 
orthogonal type, MF are in principle template-free, al- 
though in practice one can construct a template-based 
MF estimator as we have done in this paper. 

In principle, the MF are sensitive to the weighted sum 
of the bispectrum coefficients (out to the smallest scale 
measured) lOl , so the MF would naively be expected to 
contain only a subset of the same information as the bis- 
pectrum. In reality, however, this idealized expectation is 
not borne out: the bispectrum and the MF partially com- 
plement each other, and their information is not 100% 
correlated. One reason for this is the fact that the op- 
timal bispectrum estimators [111 I12j are computation- 
ally challenging to implement for current high-precision 
CMB experiments [8] [13] , and they are anyway only op- 
timal for the case of vanishing non-Gaussianity [21 [15] . 
Moreover, the bispectrum and MF are sensitive to differ- 
ent astrophysical and analysis-related systematics, given 
that they are defined in the harmonic and real space re- 
spectively. Hence, combining the constraints obtained 
by current fast though sub-optimal bispectrum estima- 
tors with those from the MF, as we do in this paper, 
provides an alternative to improving the "optimality" of 
these estimators, and makes the combined constraints 
both stronger and more robust. 

Hence an obvious question is how correlated are the MF 
and the bispectrum estimators, and consequently what is 
the combined constraint on NG from them. This is the 
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question that we address in this paper - we will show 
that the correlation between the two estimators, while 
nonzero, is far from maximal. Having calculated that, we 
compute the joint estimate of NG from both statistics. 

Minkowski junctionals methodology. The three 
Minkowski functionals V* (i — 0,1,2) describe morpho- 
logical properties of the hot and cold spots in the CMB 
temperature map. The morphology of the map, and thus 
the MF, are studied by specifying a temperature thresh- 
old v = (AT/T)/a in the map, where a is the rms of 
the fractional temperature fluctuation AT/T, hereafter 
simply denoted as /. Specifically, Vo is the area frac- 
tion of the regions above the temperature threshold, V\ 
is their boundary length, and V% is the geodesic curva- 
ture integrated along their boundary, which in a com- 
pact S2 space is related to the Euler characteristic x by 
X = V2 + Vo/2-k [TBJ. The MF can be expressed as inte- 
grals of functions of the anisotropy field and its deriva- 
tives over the compact space of the CMB sky. For explicit 
expressions see e.g. jTOl |TS] ; we shall adopt these opera- 
tionally convenient forms to calculate the MF for a given 
map. 

If the temperature fluctuations are Gaussian, the en- 
semble averages of the Minkowski functionals have an- 
alytic expressions that are completely specified by the 
two-point statistics (variance) of the fluctuations, <Tq and 
a\{= (|V/| 2 )) 17 . On the other hand, when the fluctu- 
ations are weakly non-Gaussian and the cumulants (f n ) c 
(where "c" stands for the connected part) satisfy the hi- 
erarchical ordering (f n ) c ~ Oo™ -2 , one can obtain an 
order-by-order expansion in powers of o~q for the aver- 
age of the Minkowski functionals [HI US]. In this let- 
ter, we consider the first order in the hierarchical non- 
Gaussian expansion, which in addition to 00 and o~\ de- 
pends on the three-point statistics (skewness) of the field: 
S = (/ 3 ), Si ee (/ 2 V 2 /), S n ee 2(|V/| 2 V 2 /). The two 
variance and the three skewness parameters can be cal- 
culated from theory by integrating over the power spec- 
trum and bispectrum of the CMB field, respectively; for 
explicit expressions, see [HI US] ■ In the special case of the 
local- type primordial non-Gaussianity, S, Si, Sn are all 
linearly proportional to /nl- 

In this Letter, we use the co-added V+W band data 
from the WMAP seven-year results [20 to obtain our 
constraints on /nl- The V and W bands are chosen 
for they are the most foreground-free. For this purpose, 
we generate 1000 simulations of the WMAP data fol- 
lowing the procedure given in Appendix A of |21j . The 
only difference (aside from using the WMAP 7 cosmolog- 
ical model) is that we used a uniform weighting for the 
maps, rather than the slightly more complicated weight- 
ing given there, since it only gives a marginal improve- 
ment in estimating /nl- Each of our simulated map is 
the sum of three components: 1) the Gaussian CMB re- 
alizations (the "signal") based on the CMB power spec- 
trum calculated assuming the best-fit WMAP seven-year 
cosmology including the effect of beam smearing, 2) in- 
strumental noise modeled as the Poisson process with the 



rms noise per pixel a / \J A Q b s , where a is the rms noise per 
observation and -/V b s is the number of observations per 
pixel, and 3) unresolved point sources modeled as the 
Poisson realizations from assuming a single population 
of sources with a fixed frequency-independent flux whose 
flux strength and number density roughly reproduce the 
source power spectrum and bispectrum measured from 
the WMAP Q band. The latter two components are 
modeled to closely match the systematics expected in the 
V+W co-added map. We then mask both the WMAP 
data and our simulated maps by using the KQ75 mask. 

To make predictions for the ensemble average of the 
Minkowski functionals when various observational effects 
are present, we should also include these effects in the 
calculations of the two variance parameters and three 
skewness parameters. Each of these parameters has con- 
tributions from the noise part - instrumental noise and 
point sources, in addition to the beam-smeared CMB sig- 
nal part. The noise and signal contributions add up di- 
rectly since the CMB signal and noise are uncorrelated. 
We estimate these noise contributions from our simula- 
tions: we calculate the variance and skewness parameters 
for each simulated map, and take their average over the 
1000 samples; we then subtract off the signal contribu- 
tions which are known to us for these Gaussian CMB 
simulations. 

Before we proceed to the fitting procedure and obtain 
our Minkowski functional constraints on /nl, we address 
the "residual problem" in our numerical evaluation of 
the Minkowski functionals. Previous work [35] found 
that, even for a set of Gaussian CMB simulations with- 
out noise, the averages of the MF calculated for each 
map are different from their values expected from the- 
ory. As shown in [23], these residuals are generated by 
the discrete binning of the MF in the threshold v, and for 
weakly non-Gaussian maps can be calculated analytically 
and then subtracted order by order in <j . In this work, 
we instead follow [22] and calculate the residuals from 
our simulations as the difference of the sample-averaged 
means of the MF and their theoretically expected means. 
These residuals are then subtracted from the measured 
Minkowski functionals. We use the same residuals to ac- 
count for those for the non-Gaussian case: for a weakly 
non-Gaussian field, the differences are at the order of <7o- 

Before we calculate the Minkowski functionals for each 
map, we smooth the map at several different angular 
scales. This allows us to extract additional information 
from the map and tighten the constraints on /nl- Specif- 
ically, we use a Gaussian window function, and smooth 
each map at three different scales with the Full Width 
Half Maximum (FWHM) 9 set at 6 = 5', 20', 80'. Pix- 
els near the boundary of the KQ75 mask are removed 
to avoid contamination from the masked regions that 
may be introduced due to the smoothing. Then for each 
smoothed map, we calculate its three Minkowski func- 
tionals at 15 temperature thresholds from v = —3.5 to 
3.5 with equal bin size of Ai> = 0.5. 

To obtain the constraints on /nl, we perform a x 2 
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on 


/nl 




5 


-33±110 


Minkowski 


20 


-27± 52 


Functionals 


80 


-8 ± 70 




all 


-24 ± 36 


bispectrum 




46 ± 35 


MF + bisp 




12 ± 29 



TABLE I. Constraints on /nl from the CMB Minkowski func- 
tionals, bispectrum, and their combination. The analyses use 
the WMAP 7-year V+W co-added map. 6 is the angular 
scale used to smooth the map for the Minkowski functional 
analysis. 

analysis, which compares theoretical predictions at a 
given /nl to the measurements, and is calculated as 

X 2 = E [^° bS - ^(/NL)] [V° hs - ^\/nl)] , 

where i and j run over all combinations of the 15 thresh- 
olds, three orders, and three smoothing scales for the 
measured Minkowski functionals. Here V° hs are the 
"observed" numerically evaluated Minkowski functionals, 
V^ h are the theoretically expected averages for the MF 
(which are functions of /nl), and C is the covariance ma- 
trix for (Vi, Vj) which we calculate from our simulations 
as 

c ii = {(v i -{v i ))(y j -{v j ))) eim , (2) 

where the angular brackets denote averaging over the 
1000 simulated maps. We then obtain our best-fit value 
of /nl, henceforth /nl > by minimizing the \ 2 ■ 

To check that our estimator for /^f is unbiased, we 
first apply it to the 1000 simulated maps either for the 
MF measurements at each smoothing scale or their com- 
bined results. We find that the average of the best-fit 
values accurately reproduces the theoretical input in our 
simulation, i.e., /nl = 0. Next, we test our estimator on 
publicly available non- Gaussian CMB maps generated 
with the local- type NG [24] , and we again find negligible 
bias (1% or less of the true /nl) m our estimator. 

Finally, we apply our estimator on the co-added V+W 
band data from the WMAP. In Table [TJ we show the con- 
straints on /nl from smoothing the map at each of the 
three angular scales, and the joint constraints from all 
scales combined, which we quote as our final MF con- 
straints: = -24 ±36. 

Bispectrum methodology. With the MF estimator of 
/nl^ obtained, we next develop /nl P ~~ the estimator from 
bispectrum. The observed CMB bispectrum is given by 




where the matrix is the Wigner-3j symbol, and ai m is 
the spherical harmonic transform of the temperature 
anisotropy map. In the local-type NG model, B^^^ 
is linearly proportional to /nl- 

We follow the prescription that uses the KSW [35] es- 
timator to calculate /nl P from CMB maps (see also [26] 
for the exact implementation that we use). In brief, the 
KSW is a cubic (in the temperature field) estimator of 
non-Gaussianity; it is nearly minimum-variance and com- 
putationally fast, and can straightforwardly deal with 
partial sky coverage and inhomogeneous noise. The first 
ingredient in using KSW is to calculate the Fisher ma- 
trix F corresponding to /nl! for this we need the theo- 
retical bispectrum B^°^ which can be calculated with 
the help of transfer functions from CAMB [57j . Further- 
more, KSW requires filtered maps A{ft,r) and B{ft,r) 
from which the skewness S of the field can be calculated; 
these filtered maps can be computed using HEALPix (by 
way of HealPy) to perform the forwards and backwards 
spherical harmonic transforms that are necessary in their 
computation. Given the skewness and the Fisher matrix, 
the KSW estimator for /nl is 



To account for the masking of the CMB sky, we make 

the Substitution S S cut = S/f sky + Senear [2E]- ^linear 

is an addition to skewness and is calibrated to account 
for partial-sky observations 

linear / T % d,T j <fft [A {ft, r) (B^ m (ft, r)) MC 

/sky J J 

+2B{ft,r)(A sim {ft,r)B sim {ft,r)) MC } . (5) 

The subscripted filtered maps A s { m and i? s im are cre- 
ated from Python-produced Gaussian Monte Carlo real- 
izations of the cut CMB sky; the brackets ()mc indicate 
an average over 300 of the maps. The simulated maps 
were produced as outlined earlier when we discussed the 
MF. 

Applying the bispectrum/KSW estimator to the co- 
added V+W band data of the WMAP, we obtain the 
constraint on the local NG to be /nl P = 46 ± 35. 

Combined analysis. In addition to obtaining the con- 
straints on /nl separately from the MF and bispectrum 
analyses, we would like to combine them to extract a 
more stringent and robust result. To make the problem 
tractable, we opt to consistently combine the estimators 
of /nl from these two analyses, rather than attempting 
to find the covariance between the observables, i.e. the 
MF and bispectrum themselves. It is a reasonably good 
assumption that the two estimators of /nl satisfy a bi- 
variate Gaussian distribution, especially near the peak 
of the distribution (see Figure [l] below) . Let us organize 

the two estimators into a row- vector fNL = [fm!) /nl P ], 
and let C be the 2x2 covariance matrix for them. As- 
suming the underlying true value of /nl is /nl, we can 
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FIG. 1. Joint-distribution of /Jf L F and f^l P from 1000 sim- 
ulations of the WMAP data including point sources and in- 
strumental noise. We find a correlation coefficient of r = 
0.33 ± 0.03. The /^ L F estimates are obtained from the com- 
bination of smoothing the maps at 9 = 5', 20', 80'. 



write down the following joint-distribution for the two 
/ml estimators 



rocicr 1 / 2 



exp 



-^( f NL - fNL 



CT^fl 



NL 



fNL) T 



(6) 

where fNL = [/nl, /nl] = /nlI- Given a measurement 
of fNL, a best estimate for /nl can be obtained by max- 
imizing C Assuming that the covariance matrix does 
not depend on /nl, we find the following expressions for 
the best estimate and variance of /nl from the combined 
analysis 



INL 



1 



(NL 



1-1 ITT 



"!-l tr ■ 



(7) 



At the same time, by evaluating both /^f and /nl P 
for the 1000 simulated WMAP maps, we numerically ob- 
tain their joint distribution, as shown in Figure [T] From 



this distribution, we can deduce their correlation. We 
find that the two estimators of /nl are positively corre- 
lated, with a correlation coefficient of r = 0.33 ± 0.03. 
The covariances or off-diagonal elements of C are then 
C12 = C21 = r\/C\\C22\ recall that we already found the 
variances to be Cn = 36 2 and C22 = 35 2 . 

Using the numerically derived covariance matrix C, 
together with our best-fits for and / N l P , we find 

through Eq. Q the combined constraint to be 



/ NL = /M F + b -P = l2±29. 



(8) 



Conclusions. We evaluated, for the first time, the full 
covariance matrix for the Minkowski functional estima- 
tor of the local-type primordial non-Gaussianity /M F and 



bisp 
NL ■ 



/NL 

We found the correla- 



tive bispectrum estimator / 
tion coefficient r = 0.33 ± 0.03, and used it to combine 
the constraints from the MF and bispectrum (and their 
respective variances) to obtain the constraint in Eq. (JsJ) . 
One convenient feature of this work is that, by combining 
the constraints at the level of MF and bispectrum estima- 
tors, we make the problem tractable: an obvious first ap- 
proach could be to calculate the covariance between the 
observed bispectrum and MF themselves, but this is ex- 
tremely complicated, given that the MF and bispectrum 
are functions of many scales and/or thresholds. Com- 
bining the different estimators numerically, as we have 
done here for the case of local NG, can in principle be 
rather straightforwardly extended to other types of NG 
and other cosmological probes. This type of approach 
is therefore likely to become more widespread with new 
and better data. Combining the constraints obtained by 
bispectrum estimators with those from the Minkowski 
functionals, as we have done in this paper, provides an 
alternative to improving the "optimality" of these esti- 
mators and leads to combined constraints that are both 
stronger and more robust. 

Acknowledgements. We thank Licia Verde for useful 
comments. We acknowledge the use of the publicly avail- 
able CAMB [27] and HEALPix [29] packages. WF is 
supported by NASA grant NNX12AC99G. WF, AB and 
DH have been supported by the DOE and NSF at the 
University of Michigan. WF and DH thank the Aspen 
Center for Physics, which is supported by the NSF Grant 
No. 1066293, for the hospitality in the summer of 2012. 



[1] E. Komatsu, |Class. Quant. Gray. 27~T2 4010 (2010) 

[2] N. Bartolo, E. Komatsu, S. Matarrese, and A. Riotto, 

Phys. Rept. 402, 103 (2004). 
[3] M. Liguori, E. Sefusatti, J. Fergusson, and E. Shellard, 

Adv.Astron. 2010, 980523 (2010). 
[4] A. P. Yadav and B. D. Wandelt, Adv.Astron. 2010, 

565248 (2010). 

[5] E. Komatsu and D. N. Spergel, Phys. Rev. D 63, 063002 
(2001). 



[6] C. L. Bennett et al. , arXiv:1212.5225 (2 012). 

[7] A. P. Yadav a nd B. D. Wandelt, |Phys.Rev.Lett. 100, 
181301 (2008)| 

[8| K. M. Smith, L. Senatore, and M. Zaldarriaga, |JCAP 
0909, 006 (2009) 



[9] C. Hikage and T Matsubara, |Mon. Not. R. Astron. Soc. 
425, 2187 (2012)[ 
[10] C. Hikage, E. Komatsu, and T. Matsubara, Astrophys. 
J. 653, 11 (2006) 



5 



[11] 
[12] 

[13] 
[14] 

[15] 

[16] 



and 



D. Babich, Phys. Rev. D72, 043003 (2005). 
P. Creminelli, A. Nicolis, L. Senatore, M. Tegmark, 
M. Zaldarriaga, |JCAP 605, 004 (2006) | 
P. Ade et al. (Planck Collaboration), (2013). 
P. Creminelli, L. Senatore, and M. Zaldarriaga, JCAP 
0703, 019 (2007). 

M. Liguori, A. Yadav, F. K. Hansen, E. Komatsu, 
S. Matarrese, e^aL, [Prrys.Rev. D7 6, 105016 (2007)] 
J. Schmalzing and K. M Gorski, |Mon. Not. R. 
tron. Soc. 297, 355 (1998)] 



As- 



[17] H. Tom ita, |Progress of The oretical Physi cs 76, 952 
(1986)| 

T. Matsubara, l Astrophys. J. 584, 1 (2003)| 

T. Matsubara, |Phys. Rev. D 81, 083505 (2010)| 

E. Komatsu et al. , Astrophys. J. Suppl. 192, 18 (2011) 



[18] 
[19] 
[20] 
[21] 



E. Komatsu et al. , Astrophys. J. Suppl. 180, 330 (2009 



SI 

09) 



[22] C. Hik age et al. , |Mon. Not. R. Astron. Soc. 389, 1439 
(2008) 



[23] E. A. Lim and D. Simon, J. Cosmol. Astropart. Phys. 1, 



048 (2012) 

F. E isner and B. D. Wandelt, |Astrophys. J. Su pp. 184, 

264 (2009)| 

[25] E. Komatsu, D. N. Spergel, and B. D. Wandelt, Astro- 
phys.J. 634, 14 (2005). 



[26] A. Bec ker and D. Huterer, Phys.R ev.Lett. 109, 121302 

I (2012)| 



[27] A. Lewis, A. C hallinor, and A. Lasenby, Astrophys. J. 
I 538, 473720 00) , 

[28] A. P. S. Yadav, E. Komatsu, B. D. Wandelt, M. Liguori, 
F. K. Hansen, and S. Matarrese, Astrophys. J. 678, 578 



(2008) 

[29] K. Go rski, E. Hivon, A. Banday, B. W andelt, F. Hansen, 
et al. , Astrophys. J. 622, 759 (2005) 



